## bcbioRNASeq 0.2.8
## class: RangedSummarizedExperiment 
## dim: 47729 4 
## metadata(29): version level ... utilsSessionInfo
##   devtoolsSessionInfo
## assays(5): counts tpm length normalized vst
## rownames(47729): ENSMUSG00000000001 ENSMUSG00000000003 ...
##   ENSMUSG00000109577 ENSMUSG00000109578
## rowData names(7): broadClass description ... geneName
##   seqCoordSystem
## colnames(4): group1_1 group1_2 group2_1 group2_2
## colData names(21): sampleName description ... x5x3Bias xGC
## ================================================================================
## Upload Dir: "/Users/travis/build/hbc/bcbio_rnaseq_output_example"
## Upload Date: 2017-05-23
## R Load Date: 2019-03-27
## Level: "genes"
## Caller: "salmon"
## Organism: "Mus musculus"
## Interesting Groups: "group"
## AnnotationHub: "AH53222"
## Ensembl Release: "87"
## Genome Build: "GRCm38"

Overview

  • Principal Investigator:
  • Researcher:
  • Experiment:

bcbio run data was imported from /Users/travis/build/hbc/bcbio_rnaseq_output_example.

Count matrices

## Saving raw_counts.rda, normalized_counts.rda, tpm.rda to /Users/travis/build/hbc/bcbio_rnaseq_output_example/report/data
## Writing raw_counts, normalized_counts, tpm to /Users/travis/build/hbc/bcbio_rnaseq_output_example/report/results/2019-03-27/quality_control

The results are saved as gzip-compressed comma separated values (CSV). Gzip compression is natively supported on macOS and Linux-based operating systems. If you’re running Windows, we recommend installing 7-Zip. CSV files can be opened in Excel or RStudio.

  • normalized_counts.csv.gz: Use to evaluate individual genes and/or generate plots. These counts are normalized for the variation in sequencing depth across samples.
  • tpm.csv.gz: Transcripts per million, scaled by length and also suitable for plotting.
  • raw_counts.csv.gz: Only use to perform a new differential expression analysis. These counts will vary across samples due to differences in sequencing depth, and have not been normalized. Do not use this file for plotting genes.

Read metrics

Total reads

High quality RNA-seq samples ideally should have at least 10 million reads per sample.

Mapped reads

The number of mapped reads should correspond to the number of total reads.

Mapping rate

The genomic mapping rate represents the percentage of reads mapping to the reference genome. Low mapping rates are indicative of sample contamination, poor sequencing quality or other artifacts.

Number of genes detected

Gene detection saturation

We should observe a linear trend in the number of genes detected with the number of mapped reads, which indicates that the sample input was not overloaded.

Exonic mapping rate

Ideally, at least 60% of total reads should map to exons.

Intronic mapping rate

The majority of reads should map to exons and not introns.

rRNA mapping rate

Samples should have a ribosomal RNA (rRNA) contamination rate below 10%.

5’->3’ bias

RNA-seq data can have specific biases at either the 5’ or 3’ end of sequenced fragments. It is common to see a small amount of bias, especially if polyA enrichment was performed, or if there is any sample degradation. If a large amount of bias is observed here, be sure to analyze the samples with a Bioanalyzer and check the RIN scores.

Counts per gene

Generally, we expect similar count spreads for all genes between samples unless the library sizes or total RNA expression are different.

We recommend visualizing counts normalized with the Trimmed Mean of M-Values (TMM) method here (Robinson and Oshlack 2010). TMM normalization equates the overall expression levels of genes between samples under the assumption that the majority of them are not differentially expressed. Therefore, by normalizing for total RNA expression by sample, we expect the spread of the TMM-normalized counts per gene to be similar for every sample.

## 29858 non-zero genes detected
## Applying trimmed mean of M-values (TMM) normalization

TPM per biotype

Different RNA-seq processing methods can preferentially capture a subset of the RNA species from the total RNA. For example, polyA selection should select for mostly coding genes and skip a large percentage of non-polyA non-coding RNA. Here we make boxplots of the TPM for the top 12 biotypes with the most genes assigned to them for each sample.

## Warning: Factor `geneBiotype` contains implicit NA, consider using
## `forcats::fct_explicit_na`

TPM per broad biotype class

The Ensembl biotype clasifications are too specific to plot them all. Here we have grouped the biotypes into broad classes and plot boxplots of the TPM for each sample.

Count density

Generally, we expect similar count spreads for all genes between samples unless the total expressed RNA per sample is different.

## 29858 non-zero genes detected
## Applying trimmed mean of M-values (TMM) normalization

Fit modeling

Several quality metrics are first assessed to explore the fit of the model, before differential expression analysis is performed.

Dispersion

The following plot shows the dispersion by mean of normalized counts. We expect the dispersion to decrease as the mean of normalized counts increases.

## Coercing bcbioRNASeq to DESeqDataSet with DESeq2 1.22.2
## converting counts to integer mode
## estimating size factors
## estimating dispersions
## gene-wise dispersion estimates
## mean-dispersion relationship
## final dispersion estimates
## fitting model and testing

Variance stabilization

These plots show the standard deviation of normalized counts using log2(), varianceStabilizingTransform(), rlog(), and tmm() by rank(mean). Note that all counts shown are log2 scale.

## Skipping regularized log

Sample similarity analysis

Before performing similarity analysis, we transform counts to log2, which acts to minimize large differences in sequencing depth and helps normalize all samples to a similar dynamic range. For RNA-seq count data, variance increases with the mean. Logarithmic transformation of normalized count values with a small pseudocount will account for large variations seen between the highest expressing genes so that these genes won’t dominate the PCA plots. However, due to the strong noise among low count values due to Poisson, the general log2 transformation will amplify this noise, and instead, low count genes will now dominate the PCA plots. So instead, we use a regularized log (“rlog”; base 2) transformation that gives similar results for high counts as a log2 transformation but also shrinks the values of low counts towards the genes’ average across samples. We do this with the rlog() function in the DESeq2 package (Love, Huber, and Anders 2014), which we will later use for differential gene expression analysis.

Principal component analysis (PCA)

PCA (Jolliffe 2002) is a multivariate technique that allows us to summarize the systematic patterns of variations in the data. PCA takes the expression levels for genes and transforms it in principal component space, reducing each sample into one point. Thereby, we can separate samples by expression variation, and identify potential sample outliers. The PCA plot is a way to look at how samples are clustering.

## Using vst counts
## Plotting PCA using 500 genes

Covariates correlation with PCs

When multiple factors may influence the results of a given experiment, it is useful to assess which of them is responsible for the most variance as determined by PCA. We adapted the method described by Daily et al. where they integrated a method to correlate covariates with principal components values to determine the importance of each factor.

Here we are showing the correlational analysis of the rlog transformed count data’s principal components with the metadata covariates of interest. Significant correlations (FDR < 0.1) are shaded from blue (anti-correlated) to orange (correlated), with non-significant correlations shaded in gray.

## 
## running pca and calculating correlations for:
## un-scaled  data in pca;
## pve >= 5%;
## kendall cor

Hierarchical clustering

Inter-correlation analysis (ICA) is another way to look at how well samples cluster by plotting the correlation between the expression profiles of the samples.

## Using vst counts

Methods

RNA-seq counts were generated by bcbio and bcbioRNASeq using Salmon (Patro et al. 2017). Counts were imported into R using tximport (Soneson, Love, and Robinson 2016) and DESeq2 (Love, Huber, and Anders 2014). Gene annotations were obtained from Ensembl. Plots were generated by ggplot2 (Wickham 2009). Heatmaps were generated by pheatmap (Kolde 2015).

R session information

YAML params

## $results_dir
## [1] "results/2019-03-27/quality_control"
## 
## $bcb_file
##                                                                       bcb 
## "/Users/travis/build/hbc/bcbio_rnaseq_output_example/report/data/bcb.rda" 
## 
## $data_dir
## [1] "data"

_setup.R

## library(bcbioRNASeq)
## library(DESeq2)
## library(DEGreport)
## library(knitr)
## library(tidyverse)
## 
## # Set seed for reproducibility
## set.seed(1454944673L)
## 
## opts_chunk[["set"]](
##     autodep = TRUE,
##     bootstrap.show.code = FALSE,
##     cache = TRUE,
##     cache.lazy = TRUE,
##     dev = c("png", "pdf"),
##     fig.height = 10L,
##     fig.retina = 2L,
##     fig.width = 10L,
##     highlight = TRUE,
##     prompt = TRUE,
##     tidy = FALSE
## )
## 
## theme_set(
##     theme_paperwhite(
##         base_size = 14L,
##         legend_position = "right"
##     )
## )

sessioninfo::session_info()

## ─ Session info ──────────────────────────────────────────────────────────
##  setting  value                       
##  version  R version 3.5.3 (2019-03-11)
##  os       macOS High Sierra 10.13.3   
##  system   x86_64, darwin15.6.0        
##  ui       X11                         
##  language (EN)                        
##  collate  en_US.UTF-8                 
##  ctype    en_US.UTF-8                 
##  tz       GMT                         
##  date     2019-03-27                  
## 
## ─ Packages ──────────────────────────────────────────────────────────────
##  package                * version   date       lib
##  acepack                  1.4.1     2016-10-29 [1]
##  affy                     1.60.0    2018-10-30 [1]
##  affyio                   1.52.0    2018-10-30 [1]
##  annotate                 1.60.1    2019-03-07 [1]
##  AnnotationDbi            1.44.0    2018-10-30 [1]
##  AnnotationFilter         1.6.0     2018-10-30 [1]
##  AnnotationHub            2.14.2    2018-12-04 [1]
##  assertive.base           0.0-7     2016-12-30 [1]
##  assertive.code           0.0-3     2018-10-21 [1]
##  assertive.data           0.0-3     2018-11-21 [1]
##  assertive.files          0.0-2     2016-05-10 [1]
##  assertive.numbers        0.0-2     2016-05-09 [1]
##  assertive.properties     0.0-4     2016-12-30 [1]
##  assertive.sets           0.0-3     2016-12-30 [1]
##  assertive.strings        0.0-3     2016-05-10 [1]
##  assertive.types          0.0-3     2016-12-30 [1]
##  assertthat               0.2.1     2019-03-21 [1]
##  backports                1.1.3     2018-12-14 [1]
##  base64enc                0.1-3     2015-07-28 [1]
##  basejump               * 0.7.2     2019-01-27 [1]
##  bcbioBase              * 0.4.1     2019-01-27 [1]
##  bcbioRNASeq            * 0.2.8     2019-03-27 [1]
##  Biobase                * 2.42.0    2018-10-30 [1]
##  BiocGenerics           * 0.28.0    2018-10-30 [1]
##  BiocManager              1.30.4    2018-11-13 [2]
##  BiocParallel           * 1.16.6    2019-02-10 [1]
##  biomaRt                  2.38.0    2018-10-30 [1]
##  Biostrings               2.50.2    2019-01-03 [1]
##  bit                      1.1-14    2018-05-29 [1]
##  bit64                    0.9-7     2017-05-08 [1]
##  bitops                   1.0-6     2013-08-17 [1]
##  blob                     1.1.1     2018-03-25 [1]
##  broom                    0.5.1     2018-12-05 [1]
##  cellranger               1.1.0     2016-07-27 [1]
##  checkmate                1.9.1     2019-01-15 [1]
##  circlize                 0.4.5     2018-11-21 [1]
##  cli                      1.1.0     2019-03-19 [1]
##  cluster                  2.0.7-1   2018-04-13 [2]
##  codetools                0.2-16    2018-12-24 [2]
##  colorspace               1.4-1     2019-03-18 [1]
##  ComplexHeatmap           1.20.0    2018-10-30 [1]
##  ConsensusClusterPlus     1.46.0    2018-10-30 [1]
##  cowplot                  0.9.4     2019-01-08 [1]
##  crayon                   1.3.4     2017-09-16 [1]
##  curl                     3.3       2019-01-10 [1]
##  data.table               1.12.0    2019-01-13 [1]
##  DBI                      1.0.0     2018-05-02 [1]
##  DEGreport              * 1.18.1    2019-01-04 [1]
##  DelayedArray           * 0.8.0     2018-10-30 [1]
##  DESeq2                 * 1.22.2    2019-01-04 [1]
##  digest                   0.6.18    2018-10-10 [1]
##  dplyr                  * 0.8.0.1   2019-02-15 [1]
##  edgeR                    3.24.3    2019-01-02 [1]
##  ensembldb                2.6.7     2019-02-28 [1]
##  evaluate                 0.13      2019-02-12 [1]
##  forcats                * 0.4.0     2019-02-17 [1]
##  foreign                  0.8-71    2018-07-20 [2]
##  Formula                  1.2-3     2018-05-03 [1]
##  genefilter               1.64.0    2018-10-30 [1]
##  geneplotter              1.60.0    2018-10-30 [1]
##  generics                 0.0.2     2018-11-29 [1]
##  GenomeInfoDb           * 1.18.2    2019-02-12 [1]
##  GenomeInfoDbData         1.2.0     2019-03-27 [1]
##  GenomicAlignments        1.18.1    2019-01-04 [1]
##  GenomicFeatures          1.34.7    2019-03-25 [1]
##  GenomicRanges          * 1.34.0    2018-10-30 [1]
##  GetoptLong               0.1.7     2018-06-10 [1]
##  ggdendro                 0.1-20    2016-04-27 [1]
##  ggplot2                * 3.1.0     2018-10-25 [1]
##  ggrepel                  0.8.0     2018-05-09 [1]
##  GlobalOptions            0.1.0     2018-06-09 [1]
##  glue                     1.3.1     2019-03-12 [1]
##  gridExtra                2.3       2017-09-09 [1]
##  grr                      0.9.5     2016-08-26 [1]
##  gtable                   0.3.0     2019-03-25 [1]
##  haven                    2.1.0     2019-02-19 [1]
##  hexbin                 * 1.27.2    2018-01-15 [1]
##  Hmisc                    4.2-0     2019-01-26 [1]
##  hms                      0.4.2     2018-03-10 [1]
##  htmlTable                1.13.1    2019-01-07 [1]
##  htmltools                0.3.6     2017-04-28 [1]
##  htmlwidgets              1.3       2018-09-30 [1]
##  httpuv                   1.4.5.1   2018-12-18 [1]
##  httr                     1.4.0     2018-12-11 [1]
##  interactiveDisplayBase   1.20.0    2018-10-30 [1]
##  IRanges                * 2.16.0    2018-10-30 [1]
##  jsonlite                 1.6       2018-12-07 [1]
##  knitr                  * 1.22      2019-03-08 [1]
##  labeling                 0.3       2014-08-23 [1]
##  lasso2                   1.2-20    2018-11-27 [1]
##  later                    0.7.5     2018-09-18 [1]
##  lattice                  0.20-38   2018-11-04 [2]
##  latticeExtra             0.6-28    2016-02-09 [1]
##  lazyeval                 0.2.2     2019-03-15 [1]
##  limma                    3.38.3    2018-12-02 [1]
##  locfit                   1.5-9.1   2013-04-20 [1]
##  logging                  0.9-107   2019-02-10 [1]
##  lubridate                1.7.4     2018-04-11 [1]
##  magrittr                 1.5       2014-11-22 [1]
##  MASS                     7.3-51.1  2018-11-01 [2]
##  Matrix                   1.2-15    2018-11-01 [2]
##  Matrix.utils             0.9.7     2018-04-19 [1]
##  matrixStats            * 0.54.0    2018-07-23 [1]
##  memoise                  1.1.0     2017-04-21 [1]
##  mime                     0.6       2018-10-05 [1]
##  mnormt                   1.5-5     2016-10-15 [1]
##  modelr                   0.1.4     2019-02-18 [1]
##  munsell                  0.5.0     2018-06-12 [1]
##  nlme                     3.1-137   2018-04-07 [2]
##  nnet                     7.3-12    2016-02-02 [2]
##  Nozzle.R1                1.1-1     2013-05-15 [1]
##  pbapply                  1.3-4     2018-01-10 [1]
##  pheatmap                 1.0.12    2019-01-04 [1]
##  pillar                   1.3.1     2018-12-15 [1]
##  pkgconfig                2.0.2     2018-08-16 [1]
##  plyr                     1.8.4     2016-06-08 [1]
##  preprocessCore           1.44.0    2018-10-30 [1]
##  prettyunits              1.0.2     2015-07-13 [1]
##  progress                 1.2.0     2018-06-14 [1]
##  promises                 1.0.1     2018-04-13 [1]
##  ProtGenerics             1.14.0    2018-10-30 [1]
##  psych                    1.8.12    2019-01-12 [1]
##  purrr                  * 0.3.2     2019-03-15 [1]
##  R.methodsS3              1.7.1     2016-02-16 [1]
##  R.oo                     1.22.0    2018-04-22 [1]
##  R.utils                  2.7.0     2018-08-27 [1]
##  R6                       2.4.0     2019-02-14 [1]
##  RColorBrewer             1.1-2     2014-12-07 [1]
##  Rcpp                     1.0.1     2019-03-17 [1]
##  RCurl                    1.95-4.12 2019-03-04 [1]
##  rdrop2                   0.8.1     2017-09-29 [1]
##  readr                  * 1.3.1     2018-12-21 [1]
##  readxl                   1.3.1     2019-03-13 [1]
##  reshape                  0.8.8     2018-10-23 [1]
##  reshape2                 1.4.3     2017-12-11 [1]
##  rjson                    0.2.20    2018-06-08 [1]
##  rlang                    0.3.2     2019-03-21 [1]
##  rmarkdown              * 1.12      2019-03-14 [1]
##  rpart                    4.1-13    2018-02-23 [2]
##  Rsamtools                1.34.1    2019-01-31 [1]
##  RSQLite                  2.1.1     2018-05-06 [1]
##  rstudioapi               0.10      2019-03-19 [1]
##  rtracklayer              1.42.2    2019-03-01 [1]
##  rvest                    0.3.2     2016-06-17 [1]
##  S4Vectors              * 0.20.1    2018-11-09 [1]
##  scales                   1.0.0     2018-08-09 [1]
##  sessioninfo              1.1.1     2018-11-05 [1]
##  shape                    1.4.4     2018-02-07 [1]
##  shiny                    1.2.0     2018-11-02 [1]
##  stringi                  1.4.3     2019-03-12 [1]
##  stringr                * 1.4.0     2019-02-10 [1]
##  SummarizedExperiment   * 1.12.0    2018-10-30 [1]
##  survival                 2.43-3    2018-11-26 [2]
##  testthat               * 2.0.1     2018-10-13 [1]
##  tibble                 * 2.1.1     2019-03-16 [1]
##  tidyr                  * 0.8.3     2019-03-01 [1]
##  tidyselect               0.2.5     2018-10-11 [1]
##  tidyverse              * 1.2.1     2017-11-14 [1]
##  tximport                 1.10.1    2019-01-04 [1]
##  viridis                  0.5.1     2018-03-29 [1]
##  viridisLite              0.3.0     2018-02-01 [1]
##  vsn                      3.50.0    2018-10-30 [1]
##  withr                    2.1.2     2018-03-15 [1]
##  xfun                     0.5       2019-02-20 [1]
##  XML                      3.98-1.19 2019-03-06 [1]
##  xml2                     1.2.0     2018-01-24 [1]
##  xtable                   1.8-3     2018-08-29 [1]
##  XVector                  0.22.0    2018-10-30 [1]
##  yaml                     2.2.0     2018-07-25 [1]
##  zlibbioc                 1.28.0    2018-10-30 [1]
##  source                              
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Github (steinbaugh/basejump@685a16c)
##  Github (hbc/bcbioBase@f3b26d6)      
##  local                               
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.2)                      
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.3)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.3)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  Bioconductor                        
##  CRAN (R 3.5.3)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.2)                      
##  CRAN (R 3.5.0)                      
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
##  CRAN (R 3.5.0)                      
##  Bioconductor                        
## 
## [1] /Users/travis/R/Library
## [2] /Library/Frameworks/R.framework/Versions/3.5/Resources/library

utils::sessionInfo()

## R version 3.5.3 (2019-03-11)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: macOS High Sierra 10.13.3
## 
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
## 
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
## 
## attached base packages:
## [1] parallel  stats4    stats     graphics  grDevices utils     datasets 
## [8] methods   base     
## 
## other attached packages:
##  [1] hexbin_1.27.2               forcats_0.4.0              
##  [3] stringr_1.4.0               dplyr_0.8.0.1              
##  [5] purrr_0.3.2                 readr_1.3.1                
##  [7] tidyr_0.8.3                 tibble_2.1.1               
##  [9] ggplot2_3.1.0               tidyverse_1.2.1            
## [11] knitr_1.22                  DEGreport_1.18.1           
## [13] bcbioRNASeq_0.2.8           DESeq2_1.22.2              
## [15] bcbioBase_0.4.1             basejump_0.7.2             
## [17] SummarizedExperiment_1.12.0 DelayedArray_0.8.0         
## [19] BiocParallel_1.16.6         matrixStats_0.54.0         
## [21] Biobase_2.42.0              GenomicRanges_1.34.0       
## [23] GenomeInfoDb_1.18.2         IRanges_2.16.0             
## [25] S4Vectors_0.20.1            BiocGenerics_0.28.0        
## [27] rmarkdown_1.12              testthat_2.0.1             
## 
## loaded via a namespace (and not attached):
##   [1] R.utils_2.7.0                 tidyselect_0.2.5             
##   [3] RSQLite_2.1.1                 AnnotationDbi_1.44.0         
##   [5] htmlwidgets_1.3               grid_3.5.3                   
##   [7] munsell_0.5.0                 codetools_0.2-16             
##   [9] preprocessCore_1.44.0         withr_2.1.2                  
##  [11] colorspace_1.4-1              rstudioapi_0.10              
##  [13] assertive.base_0.0-7          rdrop2_0.8.1                 
##  [15] labeling_0.3                  lasso2_1.2-20                
##  [17] tximport_1.10.1               GenomeInfoDbData_1.2.0       
##  [19] mnormt_1.5-5                  bit64_0.9-7                  
##  [21] pheatmap_1.0.12               Matrix.utils_0.9.7           
##  [23] generics_0.0.2                xfun_0.5                     
##  [25] R6_2.4.0                      assertive.sets_0.0-3         
##  [27] locfit_1.5-9.1                AnnotationFilter_1.6.0       
##  [29] bitops_1.0-6                  reshape_0.8.8                
##  [31] assertthat_0.2.1              promises_1.0.1               
##  [33] scales_1.0.0                  nnet_7.3-12                  
##  [35] gtable_0.3.0                  affy_1.60.0                  
##  [37] ensembldb_2.6.7               rlang_0.3.2                  
##  [39] genefilter_1.64.0             GlobalOptions_0.1.0          
##  [41] splines_3.5.3                 rtracklayer_1.42.2           
##  [43] lazyeval_0.2.2                acepack_1.4.1                
##  [45] broom_0.5.1                   checkmate_1.9.1              
##  [47] BiocManager_1.30.4            yaml_2.2.0                   
##  [49] reshape2_1.4.3                modelr_0.1.4                 
##  [51] GenomicFeatures_1.34.7        backports_1.1.3              
##  [53] httpuv_1.4.5.1                Hmisc_4.2-0                  
##  [55] tools_3.5.3                   psych_1.8.12                 
##  [57] logging_0.9-107               affyio_1.52.0                
##  [59] assertive.strings_0.0-3       RColorBrewer_1.1-2           
##  [61] ggdendro_0.1-20               sessioninfo_1.1.1            
##  [63] Rcpp_1.0.1                    plyr_1.8.4                   
##  [65] base64enc_0.1-3               progress_1.2.0               
##  [67] zlibbioc_1.28.0               RCurl_1.95-4.12              
##  [69] prettyunits_1.0.2             rpart_4.1-13                 
##  [71] viridis_0.5.1                 pbapply_1.3-4                
##  [73] GetoptLong_0.1.7              cowplot_0.9.4                
##  [75] grr_0.9.5                     haven_2.1.0                  
##  [77] ggrepel_0.8.0                 cluster_2.0.7-1              
##  [79] magrittr_1.5                  data.table_1.12.0            
##  [81] assertive.data_0.0-3          circlize_0.4.5               
##  [83] ProtGenerics_1.14.0           hms_0.4.2                    
##  [85] mime_0.6                      evaluate_0.13                
##  [87] xtable_1.8-3                  XML_3.98-1.19                
##  [89] readxl_1.3.1                  gridExtra_2.3                
##  [91] shape_1.4.4                   compiler_3.5.3               
##  [93] biomaRt_2.38.0                crayon_1.3.4                 
##  [95] R.oo_1.22.0                   htmltools_0.3.6              
##  [97] later_0.7.5                   Formula_1.2-3                
##  [99] geneplotter_1.60.0            lubridate_1.7.4              
## [101] DBI_1.0.0                     assertive.files_0.0-2        
## [103] ComplexHeatmap_1.20.0         MASS_7.3-51.1                
## [105] assertive.numbers_0.0-2       Matrix_1.2-15                
## [107] cli_1.1.0                     vsn_3.50.0                   
## [109] assertive.types_0.0-3         R.methodsS3_1.7.1            
## [111] pkgconfig_2.0.2               GenomicAlignments_1.18.1     
## [113] foreign_0.8-71                xml2_1.2.0                   
## [115] annotate_1.60.1               XVector_0.22.0               
## [117] rvest_0.3.2                   digest_0.6.18                
## [119] ConsensusClusterPlus_1.46.0   assertive.code_0.0-3         
## [121] Biostrings_2.50.2             cellranger_1.1.0             
## [123] htmlTable_1.13.1              edgeR_3.24.3                 
## [125] curl_3.3                      shiny_1.2.0                  
## [127] Rsamtools_1.34.1              rjson_0.2.20                 
## [129] nlme_3.1-137                  jsonlite_1.6                 
## [131] viridisLite_0.3.0             limma_3.38.3                 
## [133] pillar_1.3.1                  lattice_0.20-38              
## [135] Nozzle.R1_1.1-1               httr_1.4.0                   
## [137] survival_2.43-3               interactiveDisplayBase_1.20.0
## [139] glue_1.3.1                    bit_1.1-14                   
## [141] assertive.properties_0.0-4    stringi_1.4.3                
## [143] blob_1.1.1                    AnnotationHub_2.14.2         
## [145] latticeExtra_0.6-28           memoise_1.1.0

Jolliffe, Ian. 2002. Principal component analysis. Wiley Online Library.

Kolde, Raivo. 2015. Pheatmap: Pretty Heatmaps. https://CRAN.R-project.org/package=pheatmap.

Love, Michael I, Wolfgang Huber, and Simon Anders. 2014. “Moderated Estimation of Fold Change and Dispersion for RNA-seq Data with DESeq2.” Genome Biol. 15 (12): 550. https://doi.org/10.1186/s13059-014-0550-8.

Patro, Rob, Geet Duggal, Michael I Love, Rafael A Irizarry, and Carl Kingsford. 2017. “Salmon Provides Fast and Bias-Aware Quantification of Transcript Expression.” Nat. Methods 14 (4): 417–19. https://doi.org/10.1038/nmeth.4197.

Robinson, Mark D, and Alicia Oshlack. 2010. “A Scaling Normalization Method for Differential Expression Analysis of RNA-seq Data.” Genome Biol. 11 (3): R25. https://doi.org/10.1186/gb-2010-11-3-r25.

Soneson, Charlotte, Michael I Love, and Mark D Robinson. 2016. “Differential Analyses for RNA-seq: Transcript-Level Estimates Improve Gene-Level Inferences.” F1000Res. 4 (December). https://doi.org/10.12688/f1000research.7563.1.

Wickham, Hadley. 2009. Ggplot2: Elegant Graphics for Data Analysis. Use R. Springer New York. https://doi.org/10.1007/978-0-387-98141-3.